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Abstract: We study a holographic superconductor model with momentum relaxation due 
to massless scalar fields linear to spatial coordinates (ifi = /36ux l ), where /3 is the strength of 
momentum relaxation. In addition to the original superconductor induced by the chemical 
potential^) at fi = 0, there exists a new type of superconductor induced by /3 even at // = 0. 
It may imply a new ‘pairing’ mechanism of particles and antiparticles interacting with /?, 
which may be interpreted as ‘impurity’. Two parameters /i and /? compete in forming a 
superconducting phase. As a result, the critical temperature behaves differently depending 
on /5//u.. It decreases when jH/fi is small and increases when /3/fj, is large, which is a novel 
feature compared to other models. After analysing ground states and phase diagrams for 
various /3//r, we study optical electric(cr), thermoelectric(a), and thermal(fc) conductivities. 
When the system undergoes a phase transition from a normal to a superconducting phase, 
1/uj pole appears in the imaginary part of the electric conductivity, implying infinite DC 
conductivity. If /3//r < 1, at small u, a two-fluid model with an imaginary 1/u pole and 
the Drude peak works for a, a, and R, but if ff/ fi > la non-Drude peak replaces the Drude 
peak. It is consistent with the coherent/incoherent metal transition in its metal phase. 
The Ferrell-Glover-Tinkham (FGT) sum rule is satisfied for all cases even when fi = 0. 
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1 Introduction 

Holographic methods (gauge/gravity duality) have provided novel tools to study diverse 
strongly correlated systems [1-4]. In particular, after a pioneering model of superconductor 
in holographic methods by Hartnoll, Herzog, and Horowitz(HHH) [5, 6], there have been 
extensive development of the model. We refer to [2, 3, 7] for reviews and references. 

The HHH model is translationally invariant. Because a translationally invariant system 
with finite charge density cannot relax momentum, the HHH model will exhibit an infinite 
electric DC conductivity even in the normal metal phase. Therefore, to construct more 
realistic superconductor models, it is important to incorporate momentum relaxation in 
the framework of holography. 

One way to include momentum relaxation is to break translational invariance by im¬ 
posing explicit inhomogeneous boundary conditions such as a spatially modulated scalar 
field or temporal U( 1) gauge field At, which mimicks an ionic lattice [8, 9]. These mod¬ 
els successfully yield a finite DC conductivity as well as interesting features in optical 
conductivity such as a Drude-like peak at small uj and some scaling laws at intermediate 
uj. We refer to [10-12] for further development. This idea (in particular, explicit optical 
lattice) was applied to the HHH model in [13] and, interestingly, many properties of the 
bismuth-based cuprates were observed. 

In this method, however, because of inhomogeneity of dynamic fields, the equations of 
motion become complicated coupled partial differential equations(PDE). It is technically 
involved and less flexible than ordinary differential equations(ODE), though conceptually 
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clear. Therefore, it will be efficient and complementary if we can analyse the system with 
ODEs. In this line a few ideas have been proposed and developed. 

Massive gravity models [14-17] introduce mass terms for some gravitons. It breaks bulk 
diffeomorphism invariance and consequently breaks translation invariance in the boundary 
field theory. Holographic Q-lattice models [18, 19] exploit a continuous global symmetry 
of the bulk theory, where, for example, the global phase of a complex scalar field breaks 
translational invariance. Models with massless scalar fields linear in spatial coordinate [20— 
23] take advantage of the shift symmetry 1 . This model is related to Q-lattice models. For 
example, a massless complex scalar with constant <p in (2.6) of [18] gives a massless axion 
linear in a spatial direction. Also there are models utilising a Bianchi VIIo symmetry to 
construct black holes dual to helical lattices [27-29]. All these models give us ODEs and 
yield a finite DC conductivity as expected. Furthermore, for a large class of models, the 
analytic DC conductivity formulas are available in terms of the data of the black hole 
horizon [30]. 

Building on this development based on ODE, it is natural to revisit the holographic su¬ 
perconductor models. The superconductor model combined to the massive gravity models 
and Q-lattice models have been studied in [31] and [32, 33] respectively. For massless scalar 
models, an anisotropic background case with one scalar field was addressed in [34, 35]. As 
in non-superconducting cases, the properties of theses ODE-based superconductor models 
qualitatively agree to the PDE-based model with ionic lattice [13]. 

In this paper, we study a holographic superconductor model based on a massless scalar 
model for isotropic background. The model consists of two parts: The HHH action [13] and 
two massless real scalar fields [20] . The HHH action is a class of Einstein-Maxwell-complex 
scalar action with negative cosmological constant. An on shell massless real scalar field(V’) 
is linear to spatial coordinate with proportionality constant /?, for example, if) = /3x. To 
have isotropic bulk fields, the identical scalar field is introduced for every spatial direction 
in field theory so there is only one parameter, /3, which controls the strength of momentum 
relaxation. 

The HHH model without massless real scalar sector has been studied extensively. 
See [2, 3] for review. It has two phases, normal metallic phase and superconducting phase. 
Without the complex scalar field, the black hole is Reissner-Nordstrom type and the system 
is in normal metal phase. With a finite complex scalar hair, the system is in superconductor 
phase. With massless real scalar fields, a normal metal phase still exist [20, 24] and its 
thermodynamic and transport coefficients were studied: the DC electric conductivity [20], 
DC thermoelectric and thermal conductivity [30], optical electric conductivity [22], optical 
electric, thermoelectric and thermal conductivities [23]. 

Having studied the metal phase of the model, we want to investigate the supercon¬ 
ducting phase. First, we will examine the condition in which a superconductor phase may 
exist at finite /3. Second, we will study the effect of /3 on superconducting phase transition 
and the properties of superconductor. 

Whis model with analytic solutions have been reported in [24] without specific applications to 
gauge/gravity duality. In holographic context a model with only one scalar field was studied for an aniso- 
topic background in [25, 26] 
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Interestingly, we find that there exists a new type of superconducting phase even when 
H = 0 at finite (3 in addition to the original superconductivity at (3 = 0 and finite fi. i.e. 
large (3 > T induces superconductivity as large /t > T does. Because two parameters n 
and [3 compete each other in forming superconducting phase the critical temperature may 
behave differently for (3/^ > 1 and for /3//t < 1. At zero temperature the f3 suppresses the 
superconductivity. At finite temperature for small /?/// the critical temperature decreases, 
while for large (3/fj,, the critical temperature increases. It is a novel feature compared to 
other models. 

Our main tool to analyse the superconducting phase is optical conductivities: electric(cr), 
thermoelectric (a), and thermal(k) conductivities. In holographic framework, the compu¬ 
tation of these conductivities are related to the classical dynamics of three coupled bulk 
field fluctuations(metric, gauge, scalar fields). By computing the on-shell quadratic action 
for these fluctuations we can read off the retarded Green’s functions relevant to three con¬ 
ductivities. Most papers deal with only electric optical conductivity. However, for better 
understanding, it will be good to have a complete set of three conductivities: a, a , and R. 
For a class of models, analytic formulas for DC conductivities are available [30, 36], but 
not for optical conductivities. In [23], a systematic numerical method to compute all three 
conductivities in a system with a constraint were developed based on [37, 38] . The method 
was applied to the normal metallic phase of our model, producing numerical conductivities 
reliably [23], and we will use the same method for the superconducting phase in this paper 2 . 

One of the main results of [23] is numerical demonstrations of coherent/incoherent 
metal transition 3 . In [42] metal without quasi-particle at strong coupling was classified by 
two classes: coherent metal with a Drude peak and incoherent metal without a Drude peak. 
At small oj, when momentum dissipation is weak (/3 < /r), all three optical conductivities 
fit well to the Drude form modified by Kq\ 


KnJ 
1 — ILOT 


+ I<0 , 


( 1 . 1 ) 


where Kq is the contribution from pair production affected by net charge density. For 
{3 <C fj,, Kq can be ignored and a modified Drude form is reduced to the standard Drude 
from. If /3 > fi optical conductivities do not fit to (1.1) and goes to the incoherent metal 
phase 4 , which agrees to [42], In superconducting phase, we find a similar result. If (3 < ji 
all three optical conductivities fit to 



LV 


K n r 

1 — iiOT 


+ Ko, 


( 1 . 2 ) 


where K s is supposed to be proportional to superfluid density. In superconducting phase, 
it turns out that K s ^ 0 for a and k, but K s = 0 for a. 

We have confirmed numerically the Ferrell-Glover-Tinkham (FGT) sum rule is satisfied 

2 For another numerical analysis on three optical conductivities we refer to [39]. 

3 It was shown in [40] for anisotropic background. See also [41]. 

4 See [43] for discussions on universal bounds for thermoelectric diffusion constants in incoherent metal. 
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in all cases we considered. 


/ dwRe[<r n (a;) - 07 (w)] = *-K s , (1.3) 

J o+ 2 

where 07 ( 07 ) is the electric conductivity at T < T C (T > T c ). 

This paper is organised as follows. In section 2, we introduce our holographic su¬ 
perconductor model (Einstein-Maxwell-complex scalar action with negative cosmological 
constant) incorporating momentum relaxation by massless real scalar fields. Background 
bulk solutions corresponding to superconducting phase and normal phase are obtained. By 
comparing on-shell actions of both solutions, we identify the phase transition temperature 
as a function of chemical potential and momentum relaxation parameter, which yields 3- 
dimensional phase diagrams. In superconducting phase, we also compute condensates as 
a function of temperature for given chemical potential and momentum relaxation param¬ 
eter. In section 3, we compute optical electric, thermoelectric, and thermal conductivities 
in superconducting phase and normal phase. In particular, in superconducting phase, 
we discuss the effect of momentum relaxation on conductivity in several aspects such as 
the appearance of infinite DC conductivity, Drude-nature of optical conductivity in small 
frequency range, two-fluid model, and Ferrell-Glover-Tinkham(FGT) sum rule. We also 
present a general numerical method to compute retarded Green’s functions when many 
fields are coupled. In section 4 we conclude. 

Note added: After this paper was completed, we became aware of [33] which has 
overlap with ours. 


2 Metal/superconductor phase transition 


We start with the original holographic superconductor model proposed by Hartnoll, Herzog, 
and Horowitz(HHH) [5] 


'S'hhh = [ d d+1 xy/^g 

JM 


R + d(d 9 ^ _ If 2 - ID4>| 2 - m 2 |4>| 2 
L 2 4 11 


Sgh = -2 / d d xV^TK , 


IdM 


( 2 . 1 ) 

( 2 . 2 ) 


where F = d A is the field strength for a 17(1) gauge field A and 4? is a complex scalar 
field. We have chosen units such that the gravitational constant 16irG = 1. The second 
action, Sgh> Gibbons-Hawking term, which is required for a well defined variational 

problem with Dirichlet boundary conditions. 7 is the determinant of the induced metric 
at the boundary, and K denotes the trace of the extrinsic curvature. To impose a 
momentum relaxation effect, we add the action of free massless scalars proposed in [ 20 ] 


S^= [ d d+l x\f^g 

JM 



^(< 9 V >/) 2 


7=1 


(2.3) 
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(2.4) 

(2.5) 

( 2 . 6 ) 


The action Shhh + S)/> yields the equations of motion for matters 5 

X7 m F mn + iq($*D N $ - ^r) = 0 , 

(D 2 — m 2 ) $ = 0 , 

V 2 i/ji = 0 , 

where the covariant derivative is defined by Dm& = (Vm - iqAu) <J>, and the Einstein’s 
equation 


Run - \qmn (r + " \ p2 ~ \ D ®\ 2 ~ ™ 2 I $ I 2 ~ \ 

2 V L 4 2 J (2.7) 

= 2 + 2 FmqFn Q + 2 (F) M ®F) N (fr* + D N §D M &*) , 

Since we are mainly interested in 2 + 1 dimensional systems, we will set d = 3 from 
now on. In order to construct a plane(x, y)-symmetric superconducting background, we 
take the following ansatz, 

fl'T^ T ‘^ 

ds 2 = -S(r)e _x(r) dt 2 + —— + - 7 r(dx 2 + dy 2 ) , 

SM ( 2 . 8 ) 

A = At(r)dt , 4> = 4>(r), ?/>/ = = j^dux 1 , 

where non-zero At(r) is introduced for a finite chemical potential or charge density and 
non-zero 4>(r) will be a hair of black hole, yielding a finite superconducting order parameter. 
A special form of i/ji is included for momentum relaxation, where (3 may be interpreted as 
a strength of impurity. 

Plugging the ansatz (2.8) into the equations of motion (2.4)-(2.7), we have four equa¬ 
tions (2.9)-(2.12). Maxwell’s equation (2.4) yields 


A'/ + 



2g 2 4 > 2 


A t = 0, 


(2.9) 


where can be taken to be real, since r component of Maxwell’s equation implies that the 
phase of is constant. The complex scalar field equation (2.5) becomes 


4>" + 





( 2 . 10 ) 


The massless real scalar equation (2.6) is satisfied by the ansatz (2.8). The tt and rr 

5 Index convention: M , N, ■ ■ ■ = 0,1, 2, r, and fi, u, ■ ■ ■ = 0,1,2, and i,j , • • • = 1, 2. 
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components of Einstein’s equations (2.7) give 


X 1 + + 


rq 


2 A?<& 2 e x 


= 0, 


, 2 e x A' 2 
dW _i-1_ 

23 


2 _l 

Sr r 2 


S L 2 


?n 2 $ 2 q 2 A?& 2 e x 
+ —- + ~ — hi - 


-P 2 

r 2 SL 2 


( 2 . 11 ) 

( 2 . 12 ) 


We will numerically solve a set of coupled equations of two second order differential 
equations (2.9-2.10) and two first order differential equations(2.11-2.12) by integrating out 
from the horizon(r/ l ), which is defined by S(fh) = 0, to infinity. It turns out that, by the 
regularity condition at the horizon, six initial conditions are determined by three initial 
values: 

${r h ), A[(r h ), xiTh), (2-13) 

with At(rh ) = 0, for a given r^, /3, rri 2 , and L. From regularity of the Euclidean on-shell 
action the Hawking temperature(T^) is given by 


= 9 '(r/Qe-axfo) 

47T 

= 7^2 ( 12 - 2 m 2 L 2 ^r h ) 2 - 2L 2 ^ - L 2 e*^\A' t (r h )) 2 ^ . 


(2.14) 

(2.15) 


Near boundary(r —> oo) 
expansion. 


the equations (2.9)-(2.12) are solved by the following series 


X(r) 
S(r) 
A t (r ) 

4>(r) 


~ x (0) + 


(<l>(i )) 2 q^ 1 )^ 2 ) 


2 r 2 

2 ($ H )) 2 

L 2 ' 2 L 2 


+ 


3r 3 


r 

~ V7T + 


S (1) 


2 L 2 + rL 2 + ’ 


A w _ A 


(i) 


+ 


r 

$(!) $( 2 ) 
- 1 - o -b 


(2-16) 

(2.17) 

(2.18) 
(2.19) 


where we considered a case with m 2 = —2/L 2 to be concrete. The mass m 2 is related to 
the conformal dimension A of an operator in the dual field theory, m 2 = ^ , and the 

scalar field falls off as r~ A . For m 2 = —2/L 2 and d = 3, A = 1 or 2 as shown in (2.19). 
The coefficients of (2.16)-(2.19) are identified with the field theory quantities as follows. 


A^^H, Af’^p, 9 (1) -- e/2, 

=0(2) ,$(2) ~ (Q(2)), or <jd 2 ) ~ ^ ^f 1 ) (qW) , 


( 2 . 20 ) 


where /j, p , e, 0^ and are chemical potential, charge density, energy density, A = i 
operator and its source, respectively. So 3^' should vanish for condensation of See 
appendix A for details, should be set to be zero to identify the Hawking temperature 
of the black hole with the temperature of boundary field theory. Equivalently, we may 
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rescale the time 


t —»• at, e x —> a 2 e x , A t -^A t /a, (2.21) 

where a = e -x(0) / 2 , to set x^ = 0- In practical computation, this rescaling method is 
much easier since we have to shoot out from horizon. Then the field theory temperature 
is computed as 

T = e x( ° )/2 T H , (2.22) 

where Th is defined in (2.15). However, analytic formulas from here are presented by 
assuming x® = 0. 

There are two scaling symmetries of the equations of motion. The first symmetry is 

r—> a\r , (t, x, y) —> ai(t, x, y) , L —y a\L , q—tq/ai, m—tm/ai, (3—> ai/3 , 

(2.23) 

and the second one is 


r^a 2 r, (t,x,y) (t,x,y)/a 2 , A^^-a 2 A^, S-ta^S, P -t a 2 P . (2.24) 


By taking a\ = 1/L and = 1/r^ we may define new scaled variables with tilde. 


r 



(t x, y) 


7 \{t,x,y), A fl = —A fl , m = mL , 
L r h 


o b 


p 


A 

rh ' 
(2.25) 


While performing numerics we work in terms of these tilde-variables. In practice, it is 
equivlalent to set L = r# = 1 in numerical computation. However, when we interpret final 
results, we need to scale them back carefully. 

Near horizon we have seven parameters, <h( r /J, A' t (rh), x( r h )j r h, P-, an d T. Taking 
advantage of three scaling symmetries (2.23)-(2.25), we may set = L = x( r /i) = 1- We 
shoot for a given m 2 and /?, dialing <f>(l), H((l). We shoot out from horizon f# = 1 
targetting = 0 or «|( 2 ) = 0 at r = oo. Thus, the solution is given by a line in a two 
dimensional configuration space of $(1) and A! t (l). 


2.1 Normal (metallic) phase 

Here we minimally summarize the properties of the metal phase to set up the stage for our 
paper, referring to [20] for more details. 

The normal phase of the system corresponds to the solution without condensate 


= 0 , 


(2.26) 
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and the analytic solution is given by 



(2.27) 


(2.28) 


(2.29) 


(2.30) 


where mo is determined by the condition 9(r/J = 0: 


m 0 = r\ 



(2.31) 


The normal phase of the system is described by a charged black brane solution with non¬ 
vanishing /?, which describes a metal state with a finite DC conductivity. The temperature 
is given by the Hawking temperature (2.14) : 



(2.32) 


There is an important property of the geometry which we will rely on when discussing 


the instability of the metal phase. In the zero temperature limit, the near horizon geometry 
of an extremal black brane becomes AdS2 x M^ 1 with the effective radius of AdS2 given 
by 


2 = L d+1 (d ~ l)/3 2 + (d- 2) V 
' 2 d{d- 1) P 2 + (d- 2)V 


(2.33) 


Notice that even at // = 0 case, the near horizon geometry remains to be AdS2 x 1 due 
to a finite (3. In this sense, /r and f3 will play a similar role as far as instability is concerned. 


2.2 Criterion for instability and quantum phase transitions 

Before performing a full numerical analysis to search superconducting phase, we first want 
to address a simpler question: When can the normal phase( (2.26)-(2.30)) be unstable by 
small scalar field perturbations, which may develop into a hairy black hole with nonzero 
$? To address this question we perform the analysis presented in [2, 44]. 

We look for an unstable mode of = 4>(r)e~ lut of the equation (2.5) in the background 
(2.26)-(2.30): 


0 = (V M — zgA M )(V M — igA M )<I> — m 2 4> 


(2.34) 



(2.35) 
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The normal phase is unstable if there is a normalisable solution with incoming boundary 
conditions at the horizon such that uj has a positive imaginary part. Because we are 
interested in determining the critical temperature it is enough to search a static(o; = 0 ) 
normalisable mode. 

Our numerical procedure is as follows. The equation (2.35) depends only on four 
dimensionless quantities: A, q, /3/g, and T / g. For fixed /3//r, A, and q, we shoot from 
horizon to boundary keeping dialing T//r. Near boundary, cj) falls off as r A ~ 3 and r _A . We 
search for the largest value of T / g for which the coefficient of r A ~ 3 vanishes. We repeat this 
procedure for different pairs of values of A, and q. Figure 1(a) shows the line of constant 
critical temperature (T//i) in A -q space when f3 / g = 1. 

In Figure 1(a) there is a special curve(solid one) representing quantum phase tran¬ 
sition at T = 0. We may understand the quantum phase transition by considering the 
Breitenlohner-Freedman(BF) bound. First, the effective mass of scalar field, which is read 
from (2.35), near horizon at zero temperature is 

2 n 2 

mi = lim lim (m 2 + q 2 g tt A ?) = m 2 - . (2.36) 

r-tr h T—>0 v 1 _l_ P 

fi 2 

Second, the near horizon geometry of an extremal black brane is given by AdS 2 x M 2 with 
the effective radius of AdS 2 given by (2.33) 


L 


2 

2 


L 2 

6 



(2.37) 


Third, recall that, real scalar held in the AdS^+i space with the radius L^+i is unstable 
with mass(M) below the BF bound 


M 2 L 


2 _ 
d+1 — 


(P 

4 


(2.38) 


From these three data (2.36)-(2.38), we conclude that the scalar held is unstable near 
horizon if 


= m ls L l 



(2.39) 


where m 2 L 2 must be greater than — | for scalar to be stable at boundary AdS 4 space. 
Alternatively, this result can be obtained from (2.35). For the scalar held near horizon of 
the extremal black hole, (2.35) becomes 


O 777 ^ 

4>"W) + -<t>'(ri) - f-<Kri) = 0 , (2.40) 

rj r\ L 

where rj = (r — r^) and m 2 ff is defined in (2.39). The equation (2.40) is an equation for a 
real scalar held in AdS 2 spacetime with the curvature radius of unity. Thus, we reach the 
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q 


(a) Transition temperatures at /3/p = 1 



q 

(b) Quantum phase transitions 


Figure 1. Phase boundaries in q and A space. The region above every curves is the normal metal 
phase. The orange solid curve in (a) and (b) are the same. 


same conclusion as (2.39) by (2.38). 

From the inequality (2.39), we can infer that qualitatively a larger mass (or A) 6 , 
smaller q, or larger f5/[i suppresses instability (superconductivity) at zero temperature. 
Quantitative phase boundaries for several /3//r are plotted in Figure 1(b), where the region 
above (below) the curve is stable or normal phase (unstable or superconducting phase). It 
is interesting that one can tune (3/fi to trigger a qunatum phase transition. For example, 
let us consider the system with q = 3 and A = 4, which is (c) in Figure 1(b). When (5/^ is 
small the system is in normal phase, but when /?///is large the system is in superconducting 
phase. The transition occurs between [3/fj, = 1 and /3//r = 2. The system at (d) in Figure 
1(b) must be always in the normal phase at zero temperature regardless of /3//t. However 
the system with q = 2, A = 2 and q = 3, A = 2 ((a) and (b) in Figure 1(b) respectively) 
must be always in the superconducting phase at zero temperature. Notice that, in this 
case, the superconducting transition occurs even at /i = 0, i.e. /3//t = oo. It also can be 
seen in Figure 2. 

2.3 Superconducting phase 

In the previous subsection we investigated the possibility of superconducting phase at finite 
temperature and zero temperature. Based on our result on quantum phase transition, 
Figure 1(b), we may anticipate which values of parameters (q , A, /3, ^ ) allow thermal phase 
transition. In this subsection we want to confirm our anticipation by constructing explicit 
superconducting background at finite temperature. 

Our numerical analysis is performed as follows. By shooting out from horizon we 
find a numerical solutions satisfying (2.9)-(2.12) and boundary condition = 0 and 

6 The mass m 2 is related to the conformal dimension A of an operator in the dual held theory, m 2 = 
A(A—3) 

1.2 ' 
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(a) A = 2, g = 3 (b)A = 2,g = 2 (c) A = 4, q = 3 



(d) plot (a) and (b) together (e) plot (c) extended to —/? 


Figure 2. Phase diagrams for point (a),(b), and (c) in Figure 1. The meshed surface is the phase 
boundary at the critical temperature. Dark region below the surface is superconducting phase while 
region above the surface is normal phase. 

consider the condensate of the operator of dimension two, (O^ 2 )}. See (2.20). We may 
choose the different boundary condition, = 0, but we will not deal with the case in 
this paper. At high temperature we obtain only one solution, which agrees to an analytic 
solution of normal state (2.26)-(2.30). At low temperature we find another solution with 
<h ^ 0(superconducting phase) in addition to a normal state solution (2.26)-(2.30). In this 
case it turns out that the superconducting solution always has a lower grand potential and 
becomes a ground state. The phase transition is continuous at a critical temperature(T c ). 
Figure 2 shows typical examples of phase diagrams for three points (a),(b), and (c) in 
Figure 1(b). The three dimensional information in Figure 2 may be summarized in a two 
dimensional plot, for example, in the plane of dimensionless quantities T/f.3 and /i//3. In 
practice, we have obtained such two dimensional plots first and rescaled them to make three 
dimensional plots, where /r, /?, and T have the same unit of energy. Three dimensional plots 
would be more convenient to represent overall features, even though all information can be 
compressed in two dimensional plots. 

Let us start with the point (a) and (b) in Figure 1(b). They are always in supercon¬ 
ducting phase for all /3 and fi at zero temperature. As temperature increases we expect that 
the system undergoes a phase transition from superconducting phase to normal phase. Our 
numerical analysis confirms it and the phase diagram is shown in Figure 2(a)(b), where the 
meshed surface is the phase boundary at the critical temperature. Dark region below the 
surface is superconducting phase while region above the surface is normal phase. Figure 
2(a)(b) focuses on the phase structure for small /3 . In Figure 2(d) we extend (3 axis of 
Figure 2(a)(b) to larger values and combine them for comparison, where Figure 2(a) is red. 
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The red mesh is above the black mesh, which means that a large q enhances superconduc¬ 
tivity, as at the zero temperature in Figure 1. However, the phase transition line coincides 
at /z = 0, because the effect of q enters in the combination of q/i as shown in (2.35). 

Notice that the superconducting phase exists even when /z = 0. In the HHH model [6] 
the phase transition is understood as a competition between /j, and T and in this case it is 
a competition between (3 and T. The ‘pairing mechanism’ of two cases must be different, 
because when (3 = 0 it would be due to a particle-particle pair, while when /z = 0 it would 
be due to a particle-anti-particle pair also interacting with f3, which may be interpreted as 
‘impurity’ [23]. In general, at finite /z and /?, two mechanisms will compete. 

This competition is reflected on the phase boundary surfaces in Figure 2. In words, the 
dependence of the critical temperature on /3/fj, is not monotonic, the critical temperature 
decreases when /3//z is small and increases when /3//z is large. In graphics, see the line 
at /z = 2 in Figure 2(d) or the lines at /z = 0 and /z = 2 in Figure 2(a). It is different 
from the previous studies. In Q-lattice model [32, 33] and single scalar model [34] the 
critical temperature decreases as momentum relaxation effect increases while in ionic lattice 
model [13] the critical temperature increases monotonically as lattice effect increase. 

In Figure 1(b), the point (c) is different from (a) and (b), in that the system at (c) at 
zero T could be in superconducting phase or normal phase depending on the value of /3//z; 
the critical value is around /3//z ~ 1.5. It is confirmed by considering the finite temperature 
transition in Figure 2(c), where the transition close to T = 0 occurs around at the line of 
/3//z « 1.5'. If we include negative values of f3, we obtain the phase diagram Figure 2(e). It 
looks similar to the superconducting dome in cuprate superconductor phase diagram when 
we identify (3 with a doping parameter. Note that (3 is a tunable continuous free parameter 
in the solution, while m and A is fixed in the action. We also see large A suppresses 
superconductivity by comparing Figure 2(a) and (c). 

In superconducting phase there is finite condensate, 4>( 2 ), which is an order parameter. 
For example, in Figure 3, we show the condensate as a function of temperature. At the 
critical temperature condensate starts forming and increases continuously as temperature 
goes down. At very small temperature our numerics becomes not reliable and we did 
not plot in that range. The plot is for A = 2, q = 3 which is the case of Figure 2(a). 
Three values of /z//3 are chosen: /z//3 = 0(blue), fj,/f3 = l(red), and n/f3 = lO(green). 
The color here matches the color of the lines in Figure 2(a). In other words, we compute 
condensate along the vertical line(temperature) standing any point on the colored-lines in 
Figure 2(a). The condensate increases as (3 increases, which agrees to an anisotropic case 
[34], When /3 —> oo the condensate approaches to the finite upper bound, while when 
/5 —>• 0 it approaches to the lower bound. 

We finish this subsection by discussing on the on-shell action of the ground state. To 
calculate a thermodynamic potential for the black hole solutions we calculate the on-shell 

'We have not yet obtained the solution at zero T. The data for the plot is numerically computed up to 
T/n ~ 0.01. It will be interesting to construct and analyse the zero T limit solutions more precisely, for 
example, following [45]. 
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Figure 3. Condensate for three values of /z//3 and A = 2, q = 3, which is the case of Figure 2(a). 
The color here matches the color of the lines in Figure 2(a). In other words, we compute condensate 
along the vertical line(temperature) standing on the colored-lines in Figure 2(a) 



Euclidean action(S' £ ') by analytically continuing to Euclidean time(r) 

t = —IT , S E = iSren , (2.41) 

where 5 ren consists of four-dimensional Stmik an d three dimensional Sbdy : 

-SVen = 'S'HHH + Sjj; + S’gh + <5ct • (2.42) 

'-V-"---' 

= Sbulk =>S'bdy 


The first three terms are defined in (2.1), (2.3) and (2.2) respectively, and the last term is 
the counter term for holographic renormalisation [6]: 


Sct = 



+ • Vipi + (771 ($* 


n M dM$ + $n M dM®*) + r ] 2 \^\ 2 / L)J , 

(2.43) 


which cancels the divergence of the bulk action. To fix <1>( 1 ) on the boundary we choose 
( 771 , 772 ) = (0, —1) while to fix we choose ( 771 , 772 ) = (1,1). n M = (0,0,0, y/SJr)) is an 
outgoing normal vector. See appendix A for more details. 

Let us first consider the Euclidean bulk action 


5, 


bulk 


— / d X\J Lbulk : 


(2.44) 


which defines Lbuik- It can be computed following [6]. The xx-component of the Einstein 
equation gives a useful relation: 


Gxx — “7~^ (Tbulk — R) + 


(2.45) 
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where is the Einstein tensor. The trace of the Einstein equation yields 


B 2 

G\ + Gy H- 2 — R T 'Cbulk = — R • (2-46) 

Thus the bulk Lagrangian is 

£buik = -G\-G r r -^. (2.47) 

In superconducting phase 

= - J d 3 x ^ dr | (- 2 r 5 (r)e _x/2 ) -/3 2 e _x/2 j . (2.48) 

After adding S^ dy the total on-shell action is finite and reads 

S E = SU + S^dy = £ /T dr J dx'idxy W = ^ W , (2.49) 

where we consider a homogeneous and equilibrium system. Thus V 2 is the volume of the 
system and the length of the Euclidean time circle is identified with inverse temperature(T). 
W is a thermodynamic potential per unit volume: 

w = s (1) - fr k - p r dr (e-^ - 1 ) - + . (M0) 

J r h 

In our case, one of T 1 ' 1 - 1 or is zero and W is simplified 

W = gW - p 2 r h -P 2 [ dr (e~*/ 2 - l) . (2.51) 

Jr h X ' 

In normal sate, where 9^ = —mo and % = 0, (2.51) becomes 

W=^(-4r 2 -/r 2 -2/3 2 ) . (2.52) 

With the relations for energy density(e), charge density (p), and entropy density(s) 

e = (T tt ) = 2 mo , p = (J*) = pr h , s = Airrl , (2.53) 

derived from (A. 10) and the definition of the Bekenstein-Hawking entropy, the equation 
(2.52) becomes the expression for the grand potential(O) 

W = e-Ts- pp = n, (2.54) 

which agrees to [ 20 ]. 


- 14 - 



In superconducting phase, 

e = (T tt ) = - 2S (1) , p = {j t )=A s = 4vrri, (2.55) 

where and are numerical values. Lack of the analytic relation such as (2.55) 
make it difficult to check if W = fl. Therefore, we numerically checked and found that 
fl = e — Ts — /ip = gW — fi 2 rh ^ W. However, as far as the phase diagram is concerned, 
this difference does not matter. It turns out W < H in superconducting phase so we may 
use H as our criteria for phase transition. To study thermodynamical quantities it will be 
important to understand physics of the difference between W and H. We leave it for future 
study. 

3 Optical conductivity 

In this section we study electric(cr), thermoelectric(o;), thermal(K) conductivity by con¬ 
sidering small fluctuations of relevant gauge, metric, scalar fields around the normal and 
superconducting background we obtained in the previous section. From here on, we set 
L = 1 and use the scaled variables (2.25) without tilde. 

3.1 Fluctuations for optical conductivity: equations and on-shell action 

Electric conductivity is related to a small bulk gauge field fluctuation 5A x (t, r ) 

/ °o J. . 

—e-^a x ( w,r), (3.1) 

-OO 

of which boundary dual operator is electric current. The fluctuation is chosen to be inde¬ 
pendent of x and y, which is allowed since all the background fields affecting the equations 
of motion are independent of x and y. Because of rotational symmetry in x-y plane, it is 
enough to consider 5A X . The gauge field fluctuation^Aj, (t, r)) sources to metric((5^ x (t, r)) 
and scalar field(<5?/ , i(L r)) fluctuation 

/ °o H 

—e~ tut r 2 h tx (u,r) , (3.2) 

-OO 

/ OO J. . 

—e-^(o;,r), (3.3) 

and all the other fluctuations can be decoupled. Notice that ht x (u, r) is defined to approach 
constant as r goes to infinity. 

In momentum space, the linearized equations of motion around the background (2.8) 
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are derived from (2.4)-(2.7): 


a" + 


s 


x_ 

2 


aL + 


u 


=x __ 


2 q 2 $ 2 

~S~ 


r 2 e x A' t , 

a x H g h tx — 0 , 


7 / -A* i/3 Se x , 

hf, r . H r a x H n-£ — 0 > 


r 2 w 


oj 2 e x 


tn,{ S' x' 2 \ , i/3we x 
£ + '~S _ ~2~ + ri^ g2“ /lte + g2 


-£ = o, 


(3.4) 

(3.5) 

(3.6) 


where SjXj^V*^ are background field obtained in the previous section. For normal phase 
we have the analytic solutions ((2.26)-(2.30)) but for superconducting phase we have it 
numerically. 

We solve these equations with two boundary conditions: incoming boundary conditions 
at the black hole horizon and the Dirichlet boundary conditions at the boundary. First, 
near the black hole horizon (r —> 1 ) the solutions are expanded as 

h tx = (r - ir +1 (hg + hg T \r - 1) + ■ ■ ■ ) , (3.7) 

a x = (r- l) u± (aW + (r - 1) + • • •) , (3.8) 

£ = (r-ir ± (£( / )+£( // )(r-l) + -..) , (3.9) 

where v± = and the incoming boundary condition corresponds to v = u + . Next, 

near the boundary (r —> oo) the asymptotic solutions read 

h - V°) 4- —h (2) 4- — /, (3) -4 
,l tx — 'Hx ' 2 r 3 n tx ' ' > 

a x = c 4 0) + * + • • • , 

£ = £ (0) + 4£ (2) + 4£ (3) + • • ‘ , 

and we fix the values of the leading terms as boundary conditions. 

Plugging the solutions into the renormalized action (2.42), we have a quadratic order 
on-shell action 


(3.10) 

(3.11) 

(3.12) 


S (2) = / d 3 x 


h tx (-^e^OxA't + + Ve x/ 2 /i( 


tx 


+ e xjl2 h 2 x ( 2 r 2 — V= + 


2 r 4 r 4( h 2 \ 


vs 2Vs; 


(3.13) 


f3£r 2 e x / 2 h tx 1 _x , l P(t 2 , £r 2 e x/2 £ 

2 VS 2 J 2^ 2 VS 


where we discarded the contribution from the horizon as the prescription for the retarded 
Green’s function [46]. In particular, with the spatially homogeneous ansatz (3.1)-(3.3), the 
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quadratic action in momentum space yields 

sgi + «) - 3 hgJhS + , (3.14) 

where V2 is the two dimensional spatial volume f dxdy and we omit the term proportional 
to <1since we are studying the case with = 0. The range of co is chosen 
to be positive following the prescription in [46]. 

The on shell action (3.14) plays a role of the generating functional for two-point Green’s 
functions sourced by and We may simply read off part of the two point 

functions from the first two terms in (3.14). The other three terms are nontrivial and 
we need to know the dependence of {a^, h^), on {ai°\ h^), £^}. However, thanks 
to linearity of equations (3.4)-(3.6), we can always find out the linear relation between 
and {a { x\h[ 0 J , £(°)}. We will first explain our numerical method to find 
such a relationship in a more general setup in the following subsection and continue the 
computation in that setup. 


3.2 Numerical method 


A systematic numerical method for a system with multi fields and constraints were devel¬ 
oped in [23] based on [37, 38]. We summarise it briefly and refer to [23] for more details. 

To develop a systematic method in a general setup, let us start with N fields 4> a (x, r), 
(a = 1,2, ,N), which satisfy a set of coupled N independent second order diffrential 

equations: 


<F a (a;,r) 


d d k 
(2ir) d 


e ~ ikx r p <$> a k { r ), 


(3.15) 


where the index a may include components of higher spin fields. For convenience, r p is 
multiplied such that the solution $^(r) goes to constant at boundary. For example, p = 2 
in (3.2). 

Near horizon(r = 1), solutions can be expanded as 


$ a (r) = (r - \) Ua± (ip a + (p a {r - 1) H-) , 


(3.16) 


where we omitted the subscript k for simplicity and u a ± correspond to incoming/outgoing 
boundary conditions. To compute the retarded Green’s function we choose the incoming 
boundary condition [46], fixing N initial conditions. The other N initial conditions, denoted 
by (z = 1,2,--- , N), can be chosen, for example, as 


<Pi V 2 ¥>3 


= 


1 

1 

1 ... 

1 \ 

1 

-1 

1 ... 

1 

1 

1 

. 1 

1— 1 

1 

V 1 

1 

1 ... 

T—1 

•• 1 


(3.17) 
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Every column vector yields a solution, denoted by 3>“(r), which is expanded as 

(r) + ''' + + • • • (near boundary), (3.18) 

where §“ are the sources( leading terms) of z-th solution and Of are the operator expectation 
values corresponding to sources (<5 a > 1). Notice that § and Q can be written as regular 
matrices of order N, where the superscript a runs for row index and the subscript i runs 
for column index. 

Since all N solutions {<3?f} becomes a basis set, a general solution yields 

<h a (r) = 4>f (r)c* —> S“c* + • • • + — j— + • • • (near boundary) (3.19) 

= J a + --- + ^ r + ---, (3.20) 

yOa 

with real constants c*’s. For any given J a we can always find c l 

c* = (S -1 )^ J a , (3.21) 


and the corresponding response R a is expressed as 

R a = 0“c* = Of (S -1 )j J b . (3.22) 


A general on-shell quadratic action in momentum space has the form of 



d d k 
(2vr) d 


J a _ k A ab (k)J b k + r k M ab (k)R b k 


(3.23) 


where A and B are regular matrices of order N. In matrix notation, can be understood 
as a row matrix. For example, the action (3.14) can be written in the form (3.23) with 


J a = 



R a = 



A = 


/0 -p 0 N 

o 2gW o 
o o, 



where the index ui is suppressed. With (3.22) the action (3.23) becomes 



1 f d d k 

2 L> o ( 2 ^ 

1 r d d k 

2 L> o W? 


r k [ A ab (k) + M ac o c i (s~ 1 )i(k)\4 

ja s'iR jb 
J -k u ab J k 5 


(3.24) 


(3.25) 


where the range of u is chosen to be positive following the prescription in [46]. See also 
[29] for a careful derivation of the gauge invariant Green’s function matrix. Notice that for 
one field case without mass term, this is the well known structure of the retarded Green’s 
function: A = 0 and G R ~ O/S. 
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In summary, to compute the retarded Green’s function. We need four square matrices 
of order N (the number of fields): A, B, §, O. A and B can be read off from the action (3.23), 
after taking care of all divergences by counter terms. To construct regular matrices, S and 
Q, we solve a set of differential equations N times with independent initial conditions. The 
retarded Green’s function is schematically A + B • 0 • § _1 . There is one subtlety in our 
procedure. 

For our equations there is one subtlety caused by a symmetry of the system. Solving 
the equations near horizon with the expansions (3.7)-(3.9) we find that only two of a x , 
and htj are free. Therefore, we cannot make a complete basis to construct a general 
J a . It is due to the gauge fixing g rx = 0. However, it turns out that there is a residual 
gauge transformation keeping g rx = 0, which is generated by the vector field of which 
non-vanishing component is = ee~ lut . So we may add a constant vector (Sg) along the 

residual gauge orbit. 


Sg = (O,l,*0/o;) T , (3.26) 

since J&^gtx = —i<jor 2 £ x and L^<p = f5^ x . Notice that §q satisfies the equations of motion 
(2.4)-(2.7), since the residual gauge transformation leaves the linearised equation of motion 
invariant. Therefore, our procedure is equivalent to formally adding a constant ‘solutions’ 
of the equations to the solution set {§“}. 

With the matrices § and Q, which is numerically computed, we may construct a 3 x 3 
matrix of the retarded Green’s function. We will focus on the 2x2 submatrix corresponding 
to a 1°) and h^) in (3.10). Since a is dual to 17(1) current J x and is dual to energy- 
momentum tensor Tf X 


G n G i 2 

G 21 G 22 


g jj 

X •Jx 

G T tx J x 


a 1 } T 

Jx -L tx 

rR 

GT rp rji 

J-txJ-tX 


(3.27) 


From the linear response theory, we have the following relation between the response func¬ 
tions and the sources: 


( W \ ( On G 12 ) / oi 0) \ 

{{T t ,)) [g 21 G 22 ) (kg)) ' 


(3.28) 


We want to relate these Green’s functions to the electric (a), thermal (n), thermoelectric 
(cc, a) conductivities dehned as 


( {Jx)\ = ( o- aT\ ( E x \ 
\(Qx)J [aTRTj \-(y x T)/Tj ’ 


(3.29) 


where Q x 
As shown 
expressed 


is the heat current, E x is an electric field and V 2 .T is a temperature gradient, 
in [2, 3, 23], by taking into account a diffeomorphism invariance, (3.29) can be 
as 


f (Jx) \ 

\(T tx )-g(J x )) 


f a aT \ ( 'iA(ai 0) + iihfj) A 
\aTRT)\ iuhfj )■ 


(3.30) 
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(a) /3 /m = 0.1, T/T c = 1.52,1, 

0.94,0.76,0.37 (dotted, red, 
orange, green, blue) 


(b) $/n= 1, T/Tf = 3.2,1, 0.89, 
0.66,0.27 (dotted, red, orange, 
green, blue) 


(c) /3 /m -»■ oo(m = 0), T/T c = 
13.2,3.5,1,0.95,0.7,0.4,0.25 
(dashed, dotted, red, orange, 
green, blue, purple) 


Figure 4. Electric conductivity(cr) for three cases /3//t = 0.1,1 and oo(or fi = 0) 



Re[<7) 


w/T 



wfT 


Re[a] 




(a) /3 /m = 0.1. Data points and fitting curves (3.32) (b) /3 /m = 1- Data points and fitting curves (3.33) 
The purple line fits well too. 


Figure 5. Near ui = 0 of Figure 4(a) and (b). T/T c = the same color as Figure 4. Dots are the 
same data in Figure 4 and solid lines are Drude-like fits. 


The comparison of (3.28) and (3.30) yields 

_ iGn i(Gufi—Gi2) \ 

i(GiiM-G2i) i(G22-G22(^=0)+M(-Gi2-G 2 i+GiiM)) ) ' ( 3 ‘ 31 ) 

UJ UJ / 

3.3 Electric/thermal/thermoelectric conductivites 

Figure 4 shows examples of electric optical conductivities(cr(u;)) for three cases of /3//i: 
/?/// = 0.1,1 and, oo(n = 0). This choice of parameters also corresponds to the cases with 
the green ((3/n = 0.1), red(/?/// = 1), and blue(/3//x = oo) lines in Figure 2(a) and 3. The 
color of curves represents temperature ratio, T/T Cl where T c is the critical temperature. 
The numerical values of temperature ratio are shown in the caption. In particular the 
dotted black curve 8 is for the temperature above T c , which is in metal phase and the red 
curve corresponds to the critical temperature (in practice, it is slightly higher than the 

8 There is also a dashed grey curve in (c) at fj, = 0. It is not distinguishable from the red and dotted 
black curves in Figure 4 and 7, but distinguishable in Figure 8. 


/ a aT \ 
1 aT RT J 
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transition temperature.). The first row shows the real part of electric conductivity (Refer]) 
and the second row shows the imaginary part of electric conductivity (Im[cr]). 

One common important feature in Figure 4 is the appearance of 1 /u pole in Im[cx] 
below the critical temperature, while the disappearance of \/oj pole above the critical 
temperature. By the Kramers-Kronig relation, 1/oj pole in Im[cr] implies the delta function 
at oj = 0 in Refer], Therefore, in normal phase the DC conductivity is finite due to 
momentum relaxation and in superconducting phase the DC conductivity is infinite, which 
is one of the hallmarks of superconductor. 

Roughly at T/T c > 0.5, in addition to the delta function at oj = 0, there is still a finite 
value of DC Refer] in superconducting phase. It may be interpreted as a contribution from 
normal components in superconducting phase, implying a two-fluid model. For small oj 
there is a Drude-like peak in some cases in (a) and (b) of Figure 4. For smaller (3/p ot 
at higher temperature, the peak becomes sharper. For the sake of comparison we used a 
similar scales in (a),(b) and (c) of Figure 4, which hides the structure of (a) near oj = 0. 
Therefore we zoom in Figure 4(a) in Figure 5(a). The data points well fit to the formula 
(solid lines) 


a(uj) 


= .K*_ K n r 

oj 1 — iojt ’ 


(3.32) 


where oj = oj/ p and K s and K n are supposed to be proportional to the superfluid density 
and normal fluid density. For (3/p = 1 the formula (3.32) does not work and the data(red, 
orange, green) better fit to 


ct{oj) 



oj 


K n r 


1 — ioiT 


+ K 0 , 


(3.33) 


which is shown in Figure 5(b). Kq is related to pair creation and it was necessary also in 
metal phase. The existence of Kq is most apparent in Figure 4(c), where p = K n = 0. 
Indeed (3.32) is understood as an approximation of (3.33) when Kq is negligible compared 
to K n T. As temperature is lowered I\ n and Kq is reduced while K s is enhanced. Kq 
becomes zero at low temperature (green and blue line in Figure 4(c)), but it is possible 
that K n is finite even at zero T [13, 31, 32]. 

If (3 p 7 it was shown that the coherent metal phase becomes incoherent, where the 
Drude peak becomes a non-Drude peak [23]. The Figure 5 in [23] suggests that the critical 
(3/p is around 1/2 in metal phase. Therefore, it is suspected that, If (3 > (1/2 )p, (3.33) 
does not work in superconductor phase either. Indeed we see this is the case. When (3 = p: 
the fit of Figure 5 (b) is not as good as (a) and starts deviating from (3.33) 9 . 

As temperature goes down (T < T c ) the spectral weight of Re[cr] is reduced while K s 
of Im[<r] is enhanced. This transfer of the spectral weight to K s may be quantified by the 
Ferrell-Glover-Tinkham (FGT) sum rule: 


POO 

FGT = / dwRe[<7 n (u;) - cr s (u;)] - -K s = 0 , (3.34) 

J o+ 1 

u Thc Drude nature in superconducting phase will be related to K n rather than /x. For a better under¬ 
standing of the range of applicability of the Drude model, it is important to analyse K n more extensively. 
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FGT -0.002 
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~K 


Figure 6. The numerical errors of the FGT sum rule. FGT is defined in (3.34). Data for ^ = 0.1, 

= 1. and ^ = oo are plotted as square, circle, and triangle respectively, and colors indicates the 
temperatures with the same colors in Figure 4. 


Refer] 



lm[a] 



(a) /3 /m = 0.1, T/T c = 1.52,1, 

0.94,0.76,0.37 (dotted, red, 
orange, green, blue) 


Refer] 

30 



lm[a] 



(b) P/(JL = 1, T/T c = 3.2,1,0.89, 
0.66,0.27 (dotted, red, orange, 
green, blue) 


Re[a] 

0.010 


0.005 


- 0.005 


- 0.010 


lm[a] 

0.010 


0.005 


- 0.005 

- 0.010 

(c) /3 /m -► OO (m = 0), T/T c = 
13.2,3.5,1,0.95,0.7,0.4,0.25 
(dashed, dotted, red, orange, 
green, blue, purple) 


Figure 7. Thermoelectric conductivity^) for three cases /3/m = 0.1,1 and oo(or m = 0) 


where a s (a n ) is the conductivity at T < T C (T > T c ). a n can be taken for any temperature 
for T > T c since the spectral weight is constant in metal phase [23]. We computed (3.34) 
numerically for all cases in Figure 4 and showed that the FGT sum rule is satisfied up to 
ICC 3 in Figure 6. 

At m = 0 and /3 yC 0 (Figure 4 (c)) there is no net charge and no Drude peak. The 
plots are very similar to the case /3 = 0 at finite /x- For example, see Figure 6 in [2] 
or Figure 1 in [23], where the infinite DC conductivity (1 /cm pole in the imaginary part) 
is due to translation invariance with finite /x- Here, there is no /x and no translation 
invariance. So the delta function must have a different origin, which may be a new type 
of superconductivity. Interestingly, even in this case, the FGT the sum rule works. The 
deficit of spectral function may be interpreted as a deficit of ‘particle-anti-particle pairs’, 
which will condense. It may imply a new ‘pairing mechanism’ of particles and anti particles 
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Re[7c]//7 



Im [7 c]//7 



(a) /3 /11 = 0.1, T/T c = 1.52,1, 
0.94,0.76,0.37 (dotted, red, 
orange, green, blue) 


R e[K]liJ 



Im [k]I/j 
150 



(b) p/n=l,T/T c = 3.2,1,0.89, 
0.66,0.27 (dotted, red, orange, 
green, blue) 


RepcWJ 



lm[Jc]//3 



(c) 0/n 00 (m = 0), T/T c = 

13.2,3.5,1,0.95,0.7,0.4,0.25 
(dashed, dotted, red, orange, 
green, blue, purple) 


Figure 8. Thermal conductivity( k) for three cases /3//z = 0.1,1 and oo(or fj, = 0) 



(a) a, /3/p. = 0.1, Figure 7(a) (b) k, j3 /h = 0.1, Figure 8(a) 


Figure 9. Near ui = 0. Dots are the same data in Figure 7 and 8 and solid lines are (3.32). Colors 
represent the same temperature. 


Re[i<]//3 



Figure 10. Re[«]//3 near w = 0 at /i = 0, Figure 8 (c) 


interacting with f3, which may be interpreted as kind of ‘impurity’ [23]. 

Figure 7 and 8 show the thermoelectric conductivity (a) and the thermal conductivity 
(k) for the same parameters as Figure 4. At finite fj, ((a) and (b)) the qualitative structure 
of Re [a] and Re[«] at small cj are similar to Refer], For the sake of comparison we used a 
similar scales in (a),(b). To see the structure of (a) near u = 0 we zoom in Figure 7(a) 
and 8(a) in Figure 9. The data points fit to (3.32). At fj, = 0, the DC values of Re[k] 
decreases quickly as temperature goes down as shown in Figure 10. The thermoelectric 
conductivity a vanishes, if fi = 0, which is due to particle-hole symmetry. As u —> oo, a 
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and k approaches to 


^ 2 + / 3 2 


(3.35) 


a 


a 

T ’ 


T 


which agrees to the Ward identity [47]. 

There is no 1/w pole in Im[a] for all cases. For R, at fj, = 0 there is no 1/uj pole in Im[«] 
while in (a) and (b) there is 1/u pole only in the superconducting phase. The appearance 
of this l/uo pole in R may be understood from (3.31) as follows. In superconducting phase 
there is a pole in Im[cr] so Re[Gn(0)] ^ 0. From Figure 7 we see that Re[Gi2(0)] = 
Re[Gii(0)]/z and Re[Gi2(0)] is finite if // is finite. Then, lm[k;(0)] = ^ t 1 ^ . 

Therefore, if we subtract ^ Rc( R 2 (°) f rom Im[/c], 1/w pole is expected to be disappeared, 
which we have confirmed numerically. 


4 Conclusions 

In this paper, we studied a simple holographic superconductor model incorporating momen¬ 
tum relaxation. The model consists of two parts: The HHH model [5, 6] and massless scalar 
fields sector for momentum relaxation [20], where the strength of momentum relaxation is 
parameterised by (3. 

One of the interesting features of the model is that the existence of a new type of 
superconductor induced by f3 at [A = 0. While in the HHH model the superconducting 
phase transition is understood as a competition between /r and T, in this new case, it is 
a competition between /3 and T. The ‘pairing mechanism’ of two cases must be different. 
In the new type(/3 ^ 0 and /r = 0), there is no net charge and the mechanism will be due 
to particle-anti-particle pairs also interacting with /?, which may be interpreted as kind of 
‘impurity’ [23]. The electric optical conductivity of this new superconductor satisfies the 
FGT sum rule too. The deficit of the spectral weight may be interpreted as a deficit of 
particle-anti-particle pairs which are condensed. 

With finite /r and (3 together, two superconducting mechanisms will compete. As a 
result, the dependence of the critical temperature on (3/fi is not monotonic: the critical 
temperature decreases when f3 / /j, is small and increases when /?/yu is large. It is differ¬ 
ent from the previous studies. As momentum relaxation effect increases, in a Q-lattice 
model [32, 33] and a single scalar model [34] the critical temperature decreases while in the 
ionic lattice model [13] the critical temperature increases. The condensate has the upper 
bound when /3//r —> oo and the lower bound when /3/yn —> 0. 

We studied optical electric(er), thermoelectric(a), and themial(k) conductivities. For 
all three conductivities, at small a;, a two-fluid model with a modified Drude peak works 
if (3/[i < 1: 

• K s K n r t*-i\ 

(7, Ck, K ~ l-h --:-hA 0 , (4.1) 

U) 1 — VjJT 

where I\ s and K n are supposed to be proportional to the superfluid density and normal 
fluid density and Kq is related to pair creation. For (3/n <S 1, Kq becomes negligible. The 
restriction (3/fi < 1 for (4.1) is consistent with the result in metal phase, where coherent 
metal becomes incoherent metal when /3/fj, > 1 and the Drude peak does not work [23]. 
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However, the Ferrell-Glover-Tinkham (FGT) sum rule is satisfied for all cases regardless of 

/Va 

We have not fully analyzed the parameters, K s , K n , Kq, t, their physical meanings and 
relations. The temperature dependence of the parameters are of physical importance. For 
example, K n may be related to the energy gap as studied in [13, 31-33]. The (3 dependence 
of r is relevant to the nature of dissipation. The correct identification of the superfluid 
density, proportional to I \ s , will be essential to investigate Homes’ law [48] holographically 
[49]. It will be also useful to obtain analytic formula for DC conductivities from the horizon 
data in superconducting phase as in metallic phase [30]. While the model we considered 
shows many interesting features as metal and superconductor, it also has shortcomings. 
The electric DC conductivity in normal phase is temperature independent and the insu¬ 
lator phase is lacking. It would be interesting to consider superconducting phase without 
those shortcomings. Indeed, there is a simple generalization of the model that provides a 
temperature dependent DC conductivity and an insulating phase at small temperature [50], 
so it would be interesting to construct a superconductor model based on this background. 
It would be also interesting to extend our model to the d-wave superconductors [51, 52] 
and to consider dynamical gauge fields [53]. 
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A One point functions 

We briefly summarize how to compute one point functions holographically. Let us consider 
an ADM decomposition as follows: 

ds 2 = 'y fiu dx fl dx l/ + N 2 dr 2 , (A.l) 

where, in our case, the shift vectors vanish and N is the lapse function given by 

The outward pointing normal vector is n M = (0, 0, 0, y/dJr)) and the extrinsic curvature 
is given as K /IM = 

Let us consider a renormalised action(S' ren ) consisting of a bare action(S'o) and a 
counter action(S' c t) by taking into account the holographic renormalization: 

Sren = Sq + <Sct , (A.2) 


1 

VW) 
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where 


So = £hhh + + ^gh 

-\f 2 X( a br) 2 - l 1 ^! 2 - ™ 2 I$I 2 

Z 7=1 

and 

Sct = ieb /* d3aV=7 ( - fl[71 + b *' v "” 

+T]i(<f>*n A d A <5> + <S>n A d A ( S>*) + m ( S > *$/L^ . 

To fix <f>« on the boundary we choose ( 7 / 1 , 772 ) = (0,-1) while to fix we choose 
( 771 , 772 ) = (1,1) [ 6 ]. Variations of the on-shell action with respect to fields yield 



IP" = ^E| = ^=7 (K^ - 'f'K) + , 

w= ^\ = _^~ gF ,^ 

ttI SS ren f - 55 ct 

n = yr— = -V-ffV l/>7 + y— , 
oyjj oipi 

where the variations from the counter action are 

^ =b =7 (~W U + G Ht] ~ • Vi/y 

07 ^ V ^ 2 4 

-^(277i|d>|n A ^|cl>|+r7 2 |ch| 2 /L)) , 

SSct , — rn i 


(A.5) 

(A.6) 

(A.7) 


(A.8) 


(A.9) 


The expectation values of the energy momentum tensor, the current and the scalar opera¬ 
tors in the dual field theory can be computed as 


(Tfj, u ) = lim r-LlI^, (<P) = lim r 3 —XlP , (0/) = lim r 3 —Xll 7 , 

r—> 00 yj—'y r—>• oo w—'y r—>-oo w—'y 

(0f =2 > = lim r- 1 (-^V r T - ^=7 $/£) , (o£ =1 > = lim r“ 2 (V=7 $/£) • 

' r—> 00 x ' r—>00 v 7 

(A.10) 
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